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DOES k 2 2 VARY? 

C. A. Wagner 
Geodynamics Branch 

Trajectory Analysis and Geodynamics Division 
ABSTRACT 

An attempt has been made to find a secular drift in \ 2 2 » or the phase of 
the low order and degree portion of the geogravity field. This portion may be 
associated with mass anomalies near the core-mantle boundary. From the geo- 
magnetic evidence, such anomalies might have westward drifts on the order of 
0.5 degrees/year. Tracking data on 8 synchronous satellites over a period of 
6 years were examined for residual accelerations which might be explained by a 
drift of the \ 2 2 gravity phase angle. No conclusive movement of \ 2 2 was de- 
tected. But a measured upper bound on the drift of less than 0.05 degrees/year 
is still compatible with possible slow moving irregularities in the region of the 
core-mantle boundary. 
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DOES \ 2 2 VARY? 


INTRODUCTION 

Until fairly recently, the gravitational field associated with the earth's mass 
has been assumed to be fixed. It is obvious however, that there must be small 
but important changes of the field, in time, due to a variety of causes. 

Since 1965, measurements of satellite orbit deviations have confirmed the 
existence of the solar tide induced in the earth's body. ^ ^ These observations 
have generally agreed with previous assumptions about the overall elasticity of 
the earth. In addition, recent satellite orbit studies have revealed changes in 
the earth's oblateness associated with known rotation rate changes and seasonal 
effects perhaps due to solar heating. ^ ■ W 

But so far there have been no unambiguous measurements of field changes 
due to the dynamics of the earth's own mass redistribution. The well known 
secular variations of the magnetic field (on the order of 0.5°/yr.), the almost 
certain existence of continental drift as well as the more speculative internal 
convection currents in the mantle, all imply significant changes in the 'fixed' 
gravitational field. ^ Some of these changes (if they exist) are now capable of 
being measured by satellites. 

The only known ground based observations of geogravity changes, from long 
term sea level monitoring, are inconclusive.^ The data itself (determination of 
mean sea level from tidal records) is highly variable. The estimated trends 
yield formal standard errors an order of magnitude greater than the estimations 
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themselves. However, predicted sea level trends on the basis of a 0.2°/yr. 
westward drift of only a small fraction (^20%) of the low degree and order har- 
monics are at the level of accuracy of the observations (2 to 3 cm/yr.). The 
implication is that no more than a small fraction of the low order anomalies 
are associated with the core or the core- mantle boundary. However, the in- 
terpretation of the observations is ambiguous since changes in sea level may 
be due to a number of causes not associated with gravity drift. For example, 
continental plate motions (vertical and horizontal) on the earth’s surface are 
now thought to be of the order of 0.1 to 10 cm/yr. ^ ^ 

The use of satellite tracking data over long periods of time offers the pos- 
sibility of monitoring these gravity changes unambiguously. 

THE EXPERIMENT MODEL AND DATA ANALYSIS 
The Model 

In theory the determination of gravity changes by satellite tracking could be 
made from any well and long observed orbit. In this study the record of the 
synchronous (24 hour) orbits were examined because they have almost ideal 
geometric and dynamic properties for monitoring changes in the low order field. 

Geometrically the orbits are nearly geostationary. For long periods of time 
they have remained over about the same geographic longitudes. In principal, 
such an orbit is in effect a gravity meter set up at a single longitude for long 
periods of time. The determination of the anomalous acceleration by tracking 
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spacecraft in these orbits gives repetitive ’gravity’ measurements at the same 
locations; the ideal experimental situation for monitoring long term changes. 

Dynamically, the synchronous orbits are ideal because, being stationary, 
they are extremely sensitive to the very small anomalous accelerations from 
the gravity field. On stationary orbits these accelerations build up cumulatively. 
Being amplified, they are easily detected and therefore well measured. The 
additional benefit of the high synchronous orbit is that it is sensitive to only the 
lowest degree and order gravity anomalies. These are the only anomalies 
which, because of their long wavelengths, might originate near the core- mantle 
boundary in spite of the great overburden pressures. ^ 

Considering the geogravity field developed into an infinite series of spherical 

harmonics, the geostationary orbits are especially sensitive to (or resonant with) 

the longitude dependent terms. In fact all but about 15% of the resonant acceler- 

r_i 

ation on these satellites is accounted for by the (2,2) harmonic term/ The 
experiment in this study was to try and detect a steady shift or drift of the (2,2) 
harmonic, in particular it’s phase angle (\ 2 ), examining 6 years of tracking 

data. The amplitude of the (2,2) harmonic (J 2 ) might also be considered 

variable. But it’s variability is more difficult to justify physically since it's 
change implies stress changes inside the earth. Significant steady variation in 
/ 2 2 , on the other hand, would imply no stress changes but merely lateral flow 
of anomalies in response to existing stresses (for example, at the core mantle 
boundary). 


3 



This experiment was designed with two objectives: 

1. To measure directly any long term acceleration changes on geostationary 
orbits in the same longitude region and 

2. To measure indirectly the same long term acceleration changes over 
many longitudes with respect to a drifting \ model. 

Consider the dominant resonant acceleration in longitude on a 24 hour 

... [s] 

orbit: 


where; 


X. = (A 2 2 ) sin 2 (A. - \ 2 2 ) radians/day 2 , 


2,2 


1 2rr 2 J 


2,2 



(1 + cos I) 2 /4 


( 1 ) 


a is the orbit's semimajor axis in earth radii, I is it's inclination and J 2 2 and 
\ 2 2 are related to the conventional unnormalized gravity coefficients^ 

(C- 2 , 2 and S 2 2 ) by: 


and 


^2,2 _ J 2 , 2 COS ^ ^2,2’ 


^2,2 — J 2 , 2 S * n ^ ^2,2' 


The longitude k is defined as M + a> + N - 6 e , in terms of the conventional 
Kepler elements; mean anomaly, argument of perigee, right ascension of the 
ascending node and the Greenwich hour angle. For the 24 hour circular orbit X. 
is the actual mean geographic longitude of the satellite. If we allow only X 2 2 to 
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change in time, the velocity of this change can be found from the velocity of the 
change in X. Namely, differentiating (1) with respect to time: 


d/V dX 

= ~ 2 ^2,2 [COS 2 (k - k 2 2 )] -- , 


2,2 


or 


dk 


2,2 


dt 


dk SeC 2 O - X 2,2> 

dt 2A 2i2 


( 2 ) 


Stated another way, if A k 2 2 and AX are small changes in X 2 2 and k from an 
arbitrary time when the phase angle was k 2 2 and the acceleration was k , then 
(2) implies that: 


AX 


AX 


2,2 


2 A 


2,2 


sec 2 (X - X 2 2 ) . 


(3) 


From equation (3) we can convert measurements of residual accelerations to 
measurements of phase changes. This inferred measurement of phase change 
is the basis of the experiment. 


Data Analysis 

In the first experiment (the repetitive measurement of the acceleration at 
about the same longitude), the data consists of 28 sets of mean Kepler elements 
for Intelsat 2F3 (as determined by the Communications Satellite Corporation 
between 1967 and 1971). These are shown in Table 1. The mean longitudes 
show that the satellite was within only a few degrees of 11° west during this 
time. The elements were determined by COMSAT from radar range, elevation 
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and azimuth data over periods of about a month, taken from tracking stations in 
Maine and Europe. They are separated into 3 free drift arcs by station keeping 
maneuvers performed about a year apart. Values of C 2 and S 2 were de- 
termined from the accelerated drift in each arc by fitting the X data (essentially) 
in each arc to numerical ephemerides generated by the Rapid Orbit Analysis and 
Determination program (ROAD) j- 10 ^ The accelerated drifts themselves are shown 
in Figure 1. 

From the determined values of C 2 2 and S 2 2 in each arc, best determined 
values of X were found for an orbit whose elements were the average of the 
elements in the 3 arcs. The results of this measurement are given in Table 2. 
Considering the first arc in time as the base measurement, residual accelerations 
(Arc 2 and 3 values minus the Arc 1 value) were computed for the 2nd and 3rd arcs 
of Intelsat 2F3. From these, referred phase changes were calculated from (3) 
using a value of A 2 2 (I = 0°) of 3 x 1(T 5 radians /day 2 based on a J 2 2 of 
1.8 x 10~ 5 , and a value of 2 of -15°. This data is also listed in Table 2. 
Evidently, from the overlap of the acceleration residuals themselves, there is 
no clear secular drift at this longitude over the 4 year period. 

To test the precision of the hypothesis that X 2 does have a simple linear 
drift in time, we can fit the data in Table 2 to the function: 

A\ 2 2 = Xj + X 2 t, (4) 
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where t is the time from some convenient reference. The parameter X 2 in this 
function is the drift rate X 2 2 . 

A weighted least squares solution for X 2 2 with the data in Table 2 weighted 
according to their formal standard deviations is also given in Table 2. The 
judgement above is confirmed that no clear secular trend is seen in this data. 

But the observations are limited in time span and longitude. As a result, the 
precision and reliability of this test is not sufficient to decide whether or not low 
order and degree geogravity drift exists worldwide at the same level as geo- 
magnetic drift (^0.5°/yr). 

To broaden the scope of this measurement, and to make it more precise, 12 
additional arcs of 24 hour satellites were examined with a worldwide distribution 
of longitudes and extending back to 1965. (See Table 3.) The basic data here 
(mean Kepler elements) are the same kind as in the first (direct) experiment.^ ’^ 12 ^ 
In this indirect experiment, C 2 2 and S 2 2 were best fitted (by ROAD) to the 
mean elements in all these arcs and the mean longitude residuals (observed 
minus computed values) were further analyzed for acceleration trends. The 
fixed geopotential in this determination included harmonics to (5,5) found from 
non- synchronous satellites. This analysis consisted of fitting quadratic poly- 
nomials (inthetime)tothe residuals. The AX quantities in Table 3 are the residual 
accelerations from this analysis. The AX 2 'measurements' are reductions of the 
AX's according to (3). Once again the AX 2 2 'measurements' in the 15 arcs were 
tested for a linear time drift according to the model of (4). The results of the global 
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drift measurements are shown in Figure 2 and compared to a westward drift line 
of 0.5°/yr. through this time period. 

Clearly if any low degree and order geogravity drift occurs, it is at least an 
order of magnitude slower than the equivalent wavelengths of geomagnetic drift. 

While the previous test was quite convincing, by covering more than two longi- 
tude regions the acceleration residuals were subject to systematic errors from 
poorly determined resonant terms beyond the 2nd order. Therefore, in an attempt 
to sharpen the result, a more limited data set was chosen of the best arcs in only 
the two regions 165° - 195° east, and 10° - 15° west (see Table 4). 

This final test, while conceptually cleaner than the last, gave a somewhat 
poorer estimate because the time span was slightly less and fewer arcs were 
used. In addition, the ’best' (most precisely determined) arcs were not in this 
limited set. The results; residual accelerations, residual phase shifts, and 
estimated phase drift rate in this experiment are given in Table 4 and summar- 
ized in Figure 3. 
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SUMMARY AND CONCLUSIONS 


An examination of mean elements in 15 arcs of eight 24-hour satellites 
over a 6 year period shows no significant geopotential changes. In particular, 
the data was tested for a possible drift of \ 2 , the phase angle of the 2nd order 

tesseral gravity harmonic. Over the 6 year period, the data was only compatible 
with a westward drift more than an order of magnitude slower than the westward 
drift of similar components of the geomagnetic field. However, the precision of 
the data and the short time span was not sufficient to measure rates less than 
about 0.02 “/year. 

It is possible that the earth's core drifts (or slips) westward with respect 
to the mantle at about 0.5°/yr. (as implied by the geomagnetic field drift). If so, 
then this satellite result is only compatible with an origin for the lower order and 
degree geogravity anomalies in the slower drifting portions of core-mantle boundary 
or the mantle itself. In any case there is considerable doubt that any part of the 
highly compressed lower mantle or viscous core can support even the mild shear 

stresses implied by the low order geogravity anomalies These satellite 
measurements are certainly compatible with this view also. They only definitely 
rule out an origin for these anomalies in a significantly drifting core. 
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Table 1: 

Mean Element Observations for Intelsat 2F3 






INTELSAT 2F3-1 





NO. 

TIME (MJD) 

A(E.R.) 

E 

INCL (°) 

OMEGA (°) 

NODE(°> 

MEAN (°) 

LAMBDA (° EAST) 

1 

3«6C 7 . 024 30560 

6.61 098290 

• OCO 2459 

I • 1 270 

2eb. 78 1 5 

287. 0466 

357. 8790 


y. 

3062 1 . 58 6 6CS6'2 

6.61 109304 

• OCO 1 358 

1 .2890 

313.5616 

?t? 7. 6700 

187. 1493 


3 

3064 8, 6 6666670 

6.61 099090 

.0002560 

1 .2350 

284.5675 

?8R. 1574 

271 . 4931 

348. e546 

4 

3064*. 6^ 66667Q 

6.61 079389 

• OCO 2d .3 .1 

1 . 2 19^ 

28. •'.425? 

2 6 7. c 38 4 

275. 8003 

348. 8 1 04 

5 

39684.0502778') 

6.61 0 90*99 

.0002826 

1 • 100n 

272.3*87 

288.0572 

111.5111 

?49 • 1 4 ez 

ft 

3 0 7i 9. c. : c •: ,:o 

6.611 00906 

.0004205 

1 .0 ! 60 

2( 4.8617 

2 6 7 • 4 7 2 2 

121 .691 2 

349.3375 


.30 7fll. 1 666*570 

6.6 1 071 1 1 9 

• 000570 1 

0 . 894 f' 

290. 267 4 

2 6 7. 1 4 39 

2 1 8.5057 



3.7 01 4.04 166670 

6,61 09 OH 26 

•0005517 

0.8260 

292. 8240 

289. 9«1 1 

200.9500 

250. 3999 


30645. O ' C Cn'OO 

6,61070524 

•0005676 

0.72QO 

202.6942 

292.7730 

204.4508 


i ;■ 


6 .61071218 

.00 0 49 3 3 

•:.67io 

212.0696 

292.4055 

2 12.4459 


i i 

3<66i.c : : 300 

6.61 0 71 70 7 

.0003222 

0 .6260 

20*. 21 34 

294.5369 

2 33. 31 0 7 


1 2 

3<;qc5,c-::orc:o 

6.61 050405 

.000 3249 

0.5-57.1 

2 19. 7986 

294. 1462 

246.4002 

762.3271 




INTELSAT 2F3-; 

> 




l 

4)406, 7 7 C.-t * 5 1) 

6.61111 7 i) 6 

• i ? i 3 

o. 79 8 ) 

2 c e . 9 c 4 9 

72.7948 

1 16.06! 3 

348.7097 

2 

40413, COCCfCCC 

6.61 112193 

.0 06 156 l 

*' • 8 C 4 0 

2*5.227? 

71 . 9 C 9 2 

?79.?271 

248.6229 

3 

<*C4 ;•*.< »L0 

6.61 1 03 3 97 

• : 000652 

0 , c 7 1 0 

2 2 6 . 5 1 9 9 

7a . Cl ro 

252.21 41 

248. 3314 

6 

4 C* A 0 . 0 ' * ’ 0 6 c 0 

6.61 1065^4 

• root VM 

0 .9 150 

18 1.316-' 

75,6925 

74.0376 

148. t ce7 

5 

4«'« '3, c : *<* ir c : 

6 .ol 1 26 .) *5 

•0001941 

1 

201.51 60 

7c. 66 86 . 

86 . 9993 

347. 7480 

6 

* >5 3 0 * C Is. * *'* 3 0 

6.61 1 l 77r,s 

• r 6 0 359 6 

1 • “ 6 4 * 

2 v 6 . 0 5 6 2 

75.5901 

1 09. 84 1 6 

347.4425 

7 

4 05 ? 3 • o ' : : v* c c 

6.61 1144 )o 

.000 44 2 7 

1 • l^rO 

cl' • 291 

7 5.C761 

138.1758 

346. 6627 

8 

4 05 c g. 2 « ; % : r c o 

6,6 1 1 22 1 1 6 

• 000 36 2 H 

1 • 2800 

2 20. ae 30 

75. 9267 

241 . 3623 

346. 4761 

9 

40 o ?•'■» r • c (f < co 

6,61 1 1 4 .•£ 2 o 

.0004725 

1 • 35*0 

242.2c.4«* 

76. 92 30 

1 69. 3286 

345. 8438 

1 

4* 641. ?6&6u6 70 

6.61 1 2 H 6 8 7 

.0003921 

1 .44 JO 

27?. 626^ 

6 9 • 7 C 9 0 

36.451 1 

345.5756 





INTELSAT 2F3-3 





1 

4C648.611111IC 

6.61 109900 

•00 0 28 3 > 

1.4720 

269 ,7282 

6C. 7749 

26.91 89 

345.4650 

> 

4 C *■> 7 9 . r : r OC )D 

6. 6 1 0959 30 

•0002295 

1.5170 

29C. 9030 

70.0766 

1 69.9434 

346,0141 

- 

4 C 7 60. 25 24 7 220 

6.6 1 C 98 1 20 

•0001 124 

1 .74 20 

3 3*' , 4 6 8 6 

69.3754 

3)9.4620 

3 4 7.0 49 8 

4 

4C7C 7 . cr c err cc 

6.6 1 105500 

•0000 327 

1 • 647T- 

166.8939 

70.1119 

57.4652 

347,2557 


4 •: 5 9 6 • c : i <■ o 6 o o 

6.61 CH 79 JO 

• OCO 28 1 8 

2 .0 4 70 

20 1. 6231 

7 C • 3«8D 

l 20. 4539 

347.57(6 

6 

4i c 6 s . c -n o 

6.olll 0290 

•0003444 

2.5240 

265.7846 

69.7280 

242. 1 275 

346.5583 



Table 2: Results of Repeated Measurement of the Acceleration on Intelsat 2F3 


Time Period: 39600 - 41100 MJD 
Longitude Range: 345° - 353® East 

Average Orbit Elements: a = 6.61 E.R., e = 0.004, I = 1°, A. = 348.2° East 


Acceleration and Phase Measurements 

Arc 

A 

(10 ” 5 Rad./ Day 2 ) 

Mid Arc 
Time 
(MJD) 

Longitude 

(“East) 

AA 

(10 ~ 5 Rad./Day 2 ) 

2 

(10- 2 Rad.) 

■ 

-0.015 ± 0.010 

39750 

348.2 

0.000 ± 0.010 

0.00 ± 0.17 

2 

-0.019 ± 0.006 

40500 

348.2 

-0.004 ± 0.006 

+0.07 ± 0.10 

3 

-0.020 ± 0.011 

40900 

348.2 

-0.005 ± 0.011 

+0.08 ± 0.18 


2,2 > = +0.004 ± 0.029 °/Yr. 
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Table 3: Results of Measurements of Residual Accelerations in 15 Arcs of 24- Hour Satellites 


Arc 

Time Span 
(MJD) 

Mid 

Arc 

Time 

(MJD) 

Average 

Inclination 

n 

Longitude 
Span 
(° East) 

Average 
Longitude 
(° East) 

A\ 

(10 " 5 Rad./Day 2 ) 

2 , 2 

(10 - 2 Rad.) 

Syncom 2,8 

38815-38918 

38850 

32 

65-68 

66 

-0.023 + 0.025 

-0.48 ± 0.50 

Syncom 3, 11 

39075-39263 

39150 

0 

165-172 

168 

+0.021+0.020 

-0.35 ± 0.34 

Syncom 3, 13 

39376-39531 

39450 

1 

159-161 

160 

+0.053+0.015 

-0.88 ± 0.25 

Intelsat 2F3-1 

39607-39905 

39750 

1 

349-352 

350 

+0.008+0.010 

-0.14 ± 0.17 

Intelsat 2F4-1 

40323-40608 

40450 

1 

179-194 

184 

+0.002±0.010 

-0.07 ± 0.21 

Early Bird 1 

38897-39080 

39000 

0 

324-332 

331 

+0.000+0.015 

0.00 ± 0.28 

Intelsat 2F4-2 

40617-40782 

40700 

1 

188-195 

191 

+0.000±0.005 

0.00 ± 0.14 

ATS 3-2 

40267-40337 

40300 

0 

287-288 

387 

-0.003+0.002 

-0.12 ± 0.08 

Intelsat 2F3-2 

40406-40642 

40500 

1 

346-349 

348 

+0.001±0.010 

-0.02 ± 0.17 

Skynet 1,5 

41049-41069 

41050 

1 

50 

50 

-0.015 + 0.003 

-0.39 ± 0.08 

Skynet 1, 1 

40652-40672 

40650 

2 

39-40 

39 

-0.006+0.003 

-0.32 ± 0.16 

Early Bird 2 

39096-39219 

39150 

1 

321-331 

327 

-0.006±0.004 

+0.12 ± 0.08 

Intelsat 2F3-3 

40650-41100 

40900 

1 

345-348 

347 

+0.009±0.010 

-0.18 ± 0.17 

ATS 5-2 

40742-41200 

40950 

1 

255 

255 

+0.007 + 0.003 

+0.10 ± 0.05 

Intelsat 2F4-3 

40817-41065 

40950 

2 

189-215 

194 

-0.008±0.004 

+0.22 ± 0.13 


< 2,2 y = -0.002 ± 0.016 °/Yr. 




Table 4: Results of Measurements of Residual Accelerations 
in Satellite Arcs in 2 Longitude Regions 


Arc 

A A. 

(10 - s Rad./Day 2 ) 

A ^2, 2 

(1CT 2 Rad.) 

Mid Arc Time 
(MJD) 

Average 

Longitude 

(“East) 

Syncom 3, 11 

+0.008 ± 0.020 

+0.14 ± 0.34 

39150 

168 

Intelsat 2F4-2 

-0.002 ± 0.005 

+0.06 ±0.14 

40700 

191 

Intelsat 2F4-3 

-0.007 ± 0.004 

+0.18 ± 0.13 

40950 

194 

Intelsat 2F3-2 

-0.003 ± 0.010 

+0.05 ± 0.17 

40500 

348 

Intelsat 2F3-3 

0.000 ± 0.010 

0.00 ± 0.17 

40900 

347 

Intelsat 2F4-1 

0.014 ± 0.010 

-0.30 ± 0.21 

40450 

184 

Intelsat 2F3-1 

-0.003 ± 0.010 

0.04 ± 0.17 

39750 

350 


<A 2 2 > = +0.010 ± 0.027 */Yr. 
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Figure 1. Accelerated Drifts of Intelsat 2F3 Over a 4 Year Period 
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Figure 2. Phase Shift Measurements in 15 Arcs of 24-Hour Satellites 
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Figure 3. Phase Shift Measurements For 24 Hour Satellite Arcs in 2 Longitude Regions 




